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We present STAR measurements of azimuthal anisotropy by means of the two- and four-particle 
cumulants v 2 {v 2 {2} and i>2{4}) for Au+Au and Cu+Cu collisions at center of mass energies ^/s NN = 
62.4 and 200 GeV. The difference between V2{2} 2 and «2{4} 2 is related to v 2 fluctuations (a V2 ) and 
nonfiow (62). We present an upper limit to a V2 /v 2 . Following the assumption that eccentricity 
fluctuations a £ dominate V2 fluctuations ~ we deduce the nonfiow implied for several models 
of eccentricity fluctuations that would be required for consistency with V2{2} and t>2{4}. We also 
present results on the ratio of V2 to eccentricity. 
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I. INTRODUCTION 

In non-central heavy-ion collisions, the overlap area is 
almond shaped with a long and short axis. Secondary in- 
teractions amongst the system's constituents can convert 
the initial coordinate-space anisotropy to a momentum- 
space anisotropy in the final state [TH3]- In this case, 
the spatial anisotropy decreases as the system expands 
so that any observed momentum anisotropy will be most 
sensitive to the early phase of the evolution before the 
spatial asymmetry is smoothed [1]. Ultra-relativistic nu- 
clear collisions at Brookhaven National Laboratory's Rel- 
ativistic Heavy Ion Collider (RHIC) [5] are studied in 
part to deduce whether quarks and gluons become de- 
confined during the early, high energy-density phase of 
these collisions. Since the azimuthal momentum-space 
anisotropy of particle production is sensitive to the early 
phase of the collision's evolution, observables measuring 
this anisotropy are especially interesting. The azimuth 
angle (</>) dependence of the distribution of particle mo- 
menta can be expressed in the form of a Fourier series [6] : 
dN/dcj) cx 1+ > 2v n cos n (</> — 'J), where 'J is either the 

reaction-plane angle defined by the beam axis and the 
impact parameter vectors, or the participant plane angle 
defined by the beam direction and the minor axis of the 
overlap zone [3 13- Fluctuations in the positions of nucle- 
ons within the colliding nuclei can cause deviations be- 
tween the reaction plane angle and the participant plane 
angle and the non-sphericity of the colliding nuclei may 
also enhance this effect. When energy is deposited in the 
overlap region by a finite number of collision participants, 
the energy density will necessarily possess a lumpiness as- 
sociated with statistical fluctuations. These fluctuations 
will lead to eccentricity fluctuations which can lead to v 2 
fluctuations. By definition, the eccentricity is maximum 
when calculated with respect to the participant plane. 
This plane shifts away from the reaction plane due to fluc- 
tuations. It is expected that this larger, positive definite 
eccentricity will drive the anisotropic expansion thought 
to be responsible for v 2 [7]. The eccentricity calculated 
with respect to the participant axis is called e par t and the 
eccentricity calculated with respect to the reaction plane 
is called e st( j. 

The Fourier coefficients v n can be measured and used 
to characterize the azimuthal anisotropy of particle pro- 
duction. Measurements of v 2 [E] have been taken to in- 
dicate the matter created in collisions at RHIC behaves 
like a perfect liquid with a viscosity-to-entropy ratio near 
a lower bound rj/s > 1/4-7T derived both from the uncer- 
tainty principle and string theory [TO] . This conclusion 
is primarily based on hydrodynamic model predictions 
[El [TT] . Uncertainty about the conditions at the begin- 
ning of the hydrodynamic expansion, however, leads to 
large uncertainties in the model expectations [T2] H5] . 
Since v 2 reflects the initial spatial eccentricity of the 
overlap region when two nuclei collide, fluctuations of v 2 
should depend on fluctuations in the initial eccentricity. 



Measurements of the system-size and energy dependence 
of v 2 and i>2 fluctuations are therefore useful for under- 
standing the initial conditions of the expansion phase of 
heavy-ion collisions. 

Methods used to study v 2 [H] are based on correla- 
tions either among produced particles or between pro- 
duced particles and spectator neutrons detected near 
beam rapidity ybeam- Estimates of v 2 from produced 
particles can be biased by correlations which are not 
related to the reaction or participant plane (nonflow 
S 2 = (cos(2A</>)) — (v 2 )) an d by event-by-event fluctu- 
ations of v 2 { <J V2)- Thus, an explicit measurement of 
(1*2) would require a measurement of nonflow and fluc- 
tuations. We also note that when the definition of the 
reference frame changes, from reaction plane to partici- 
pant plane for example, each of the terms v 2 , S 2 , and a V2 
can change. The experimentally observable n-particle cu- 
mulants of v 2 (labeled v 2 {2} 2 , v 2 {4} 4 , etc.) do not, how- 
ever, depend on the choice of reference frame. It has been 
shown [T5HT7] that the various analyses of v 2 based on 
produced particles are related to the second and fourth 
v 2 cumulants «2{2} and «2{4} where these are related 
to v 2 , nonflow, and fluctuations in the participant plane 
reference frame via 

v 2 {4} 2 *(v 2 f-o-t 2 (1) 

and 

v 2 {2} 2 ~v 2 {4} 2 ^d 2 + 2a 2 V2 . (2) 

These results arise because fluctuations decrease 1*2 {4} 
but increase v 2 {2} and the approximations are valid for 
a v;/( v 2) *C 1. We will discuss the effect of this ap- 
proximation later. In case the v 2 distribution is a 2D 
Gaussian in the reaction plane, the 6-particle cumulant 
v 2 {6} and higher orders will be equal to v 2 {4} and there- 
fore will not add new information. Within the accu- 
racy of the data this has been found to be the case (i.e. 
v 2 {6} f2{4}) [HI- In this approximation for the v 2 
fluctuations [T7J, v 2 {4:} is equal to the mean v 2 relative 

to the reaction plane and ^/u|{4} + a 2 2 is the mean v 2 

relative to the participant plane. We note again that a 2 2 
is not experimentally accessible without prior knowledge 
about nonflow contributions [19] . 

In this paper we present measurements of v 2 {2} and 
v 2 {4} in Au+Au and Cu+Cu collisions at ^/s NN = 200 
and 62.4 GeV. We present w 2 {2} 2 - v 2 {4} 2 « S 2 + 2a 2 2 
(called in the literature of ot ) and derive from that up- 
per limits on o~ V2 /v 2 based on several approximations. 
The upper limit assumes that v 2 fluctuations dominate 
the sum 5 2 + 2a 2 ^. This is a robust upper limit since 
larger values of o~ V2 /v 2 would require negative values of 
nonflow contrary to expectations and to measurements 
of two-particle correlations ■ We present model com- 
parisons of eccentricity fluctuations to the upper limit 
of cr V2 /v 2 . Using the same data and then alternatively 
assuming that eccentricity fluctuations drive v 2 fluctua- 
tions, we can derive the nonflow term required to satisfy 



4 



2er 2 2 for each 



the relationship v 2 {2} 2 — w 2 {4} 2 sa 5 2 
model. The 5 2 derived in this way can be compared to 
measurements of two-particle correlations [20] to check 
the validity of the models. Finally we present the ratio 
of v 2 to the initial eccentricity from the models. 

This paper does not use the method of a global fit to a 
detailed 11-parameter model of two-particle correlations 
in relative pseudorapidity and azimuth |21j . The method 
used here requires no assumptions about the shape of 
flow fluctuations or nonflow and instead considers Fourier 
harmonics of the azimuthal distributions. 

This paper is organized as follows: Section II gives 
the experimental details and cuts for the data selection. 
Section III deals with details about the Q-Cumulants 
method and the sources of systematic errors. In Section 
IV, v 2 results used in the calculation of the nonflow and 
the upper limit on v 2 fluctuations are discussed. Section 
V shows the results for the upper limit on v 2 fluctuations 
and their comparison with the eccentricity fluctuations, 
nonflow from different models and eccentricity scaling of 
v 2 for the eccentricity from different models. 



II. EXPERIMENT 

Our data sets were collected from Au+Au and Cu+Cu 
collisions at y/s NN = 62.4 and 200 GeV detected with 
the STAR detector [22] in runs IV (2004) and V (2005). 
Charged particle tracking within pseudo-rapidity < 1 
and transverse momentum px > 0.15 GeV/c was per- 
formed with the Time Projection Chamber (TPC) [25] . 
Beam-beam Counters (BBCs) and Zero-degree Calorime- 
ters (ZDCs) were used to trigger on events. We analyzed 
events from centrality interval corresponding to 0-80% 
and 0-60% of the hadronic interaction cross-section re- 
spectively for Au+Au and Cu+Cu collisions. As in pre- 
vious STAR analyses [24], we define the centrality of an 
event from the number of charged tracks in the TPC hav- 
ing pseudorapidity \r)\ < 0.5 [25]. For the v 2 analysis we 
used charged tracks with \rj\ < 1.0 and 0.15 < pr < 2.0 
GeV/c. The lower pr cut is necessitated by the accep- 
tance of the STAR detector. We varied the upper pt cut 
between 1.5 and 3.0 GeV/c to study the effect of this cut 
on the difference v 2 {2} 2 — v 2 {4} 2 . We found that v 2 {2} 
and i>2{4} increase by roughly 5% (relative) when the up- 
per pt cut is increased from 1.5 to 3.0 GeV/c, but that 
the difference between v 2 {2} 2 and «2{4} 2 changes by less 
than 1%. Only events with primary vertices within 30 cm 
of the TPC center in the beam direction were analyzed. 
The cuts used in the analysis are shown in Table |T] 



III. ANALYSIS 

We analyzed Cu+Cu and Au+Au collisions at center 
of mass energies Js NN = 62 A and 200 GeV to study the 



cut 



value 



Pt 

n 

vertex z 
vertex x,y 
fit points 

fit points/max. pts. 

dca 

trigger 



0.15 to 2.0 GeV/c 
-1.0 to 1.0 
-30.0 cm to 30.0 cm 
-1.0 cm to 1.0 cm 
>15 
>0.52 
< 3.0 cm 
Minbias 



TABLE I: Cuts used for the selection of data. Fit points are 
the number of points used to fit the TPC track, and max. 
points are the maximum possible number for that track. 



energy and system-size dependence of v 2 , nonflow and v 2 
fluctuations. From previous studies we found that it is 
not possible to use v 2 cumulants to disentangle nonflow 
effects (correlations not related to the event plane) from 
v 2 fluctuations [TO] . We have used two methods based on 
multi-particle azimuthal correlations: Q-Cumulants |26j 
for two- and four-particle cumulants to study V2{2} and 
«2{4}, and fitting the reduced flow vector q = Qj\[M dis- 
tribution to study the multi-particle v 2 . Q = e 2l< ^ 
and M is the multiplicity. The fitting of the reduced flow 
vector distribution is described in more detail in Ref . [TO] • 
The fit parameters described in that reference, «2{qfit} 



and <r 2 yn (in this paper cr 2 ot ), can be related to ^2 {2} and 
^2{4}. In Appendix [Aj we compare the q-distribution 
and Q-Cumulants results. Based on simulations, we find 
that the q-distribution method used to study v 2 by fitting 
the distribution of the magnitude of the reduced flow vec- 
tor to a function derived from the central limit theorem 
deviates more from the input values when multiplicity 
is low. For that reason, this paper presents only results 
from the Q-Cumulants method. 

The Q-Cumulants method allows us to calculate the 
cumulants without nested loops over tracks or using gen- 
erating functions [18]. For this reason it is simpler to 
perform. The cumulants calculated in this way also do 
not suffer from interference between different harmonics 
since the contributions from other harmonics are explic- 
itly removed [26]. We directly calculate the two- and 
four- particle azimuthal correlations 



(2)n|n 

(4) 



. \Qn\ 2 -M 

" ' M(M - 1) 

|Q ra | 4 + \Q2n\ 2 ~ 2Re[Q 2n Q* n Q* 

M(M — l)(M-2)(M-3) 
2(M-2)|Q„| 2 -M(A/-3) 



(3) 



(4) 



M(M — 1)(M — 2)(M — 3) ' 
where M is the number of tracks used in the analysis and 



Qn 



M 

E< 

3 



(5) 



We evaluate the terms on the right hand side of Eq. [3] 
and Eq. [4] for each event, then take the average over all 
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events. If one applies no further weighting, the two- and 
four-particle cumulant results for v n are 



within a given centrality interval. Table III in Ap- 



V n {2} 2 



(2) 
2(2) 



2 

n.n 



(4), 



n.n\n,n- 



(6) 

(7) 



It was also proposed to use weights for each event 
within a particular centrality class based on the num- 
ber of combinations of tracks for each event [26|. This 
weighting was proposed as a method to reduce the de- 
pendence of the results on multiplicity. We find however, 
that the application of number-of-combinations weights 
makes the ^2{2} and t>2{4} results more dependent on 
the width of the multiplicity bins used to define central- 
ity in our analysis. Using number-of-combination weights 
along with centrality bins defined by number of charged 
particles will lead to results that are weighted more heav- 
ily towards the higher multiplicity side of the bins and 
that effect will be stronger for four-particle correlations 
than for two-particle correlations. We also confirmed 
with simulations that without weights, the Q-Cumulant 
results for t>2{2} and i>2{4} agree better with simulation 
inputs than when weights are applied. In this paper, 
we report results without weights according to Eqs. [3] 
through [7] This method is different from that used in 
Ref. [27]. 

The systematic uncertainties on our results were esti- 
mated by evaluating our results from two different time 
periods in the run, by varying the selection criteria on 
the tracks (specifically the distance of closest approach 
of the track to the primary vertex or DC A), from the Q- 
Cumulants acceptance correction terms, and by varying 
the pt upper limit for tracks between 1.5, 2.0, and 3.0 
GeV/c. Decreasing the DCA cut and increasing the up- 
per pj- cut both increase the average p? of the analyzed 
tracks. This leads to an increase in i>2{2} and t>2{4} 
(not considered a systematic error for those data) but 
we find that the difference between i?2{2} 2 and f2{4} 2 
is nearly unchanged. This implies that the error on 
v 2 {2} 2 — «2{4} 2 due to the exact upper and lower pt 
ranges used is small. We found no difference between 
the two run periods analyzed. The acceptance correction 
applied in the analysis changes the 200 GeV Au+Au Q- 
Cumulants ^{4} results by less than 1% for all central- 
ities while the i>2{2} results change by less than 1% for 
all centralities except the 0-5% bin where they change by 
4%, and the 5-10% bin where they change by 2%. Statis- 
tical and systematic errors are shown on all results. The 
systematic errors are shown as narrow lines with wide 
caps and statistical errors are shown as thick lines with 
narrow caps. In many cases statistical errors are smaller 
than the marker size and therefore not visible. 



IV. RESULTS 

In this paper we present our results as a function of the 
average charged particle multiplicity density (dNc^/drj) 



pendix [B] provides estimates of the number of participat- 
ing nucleons AT part and {dN^/drfli 'or the centrality inter- 
vals used in this analysis. Figure[I](left) shows v 2 {2} 2 for 
200 and 62.4 GeV Au+Au collisions for charged tracks 
with 0.15 < pt < 2.0 GeV/c. The analysis is carried out 
using either all combinations of particles, independent of 
charge (CI), or using only like-sign pairs (LS). When com- 
paring the LS and CI results, we note that the LS results 
are systematically lower than the CI results for all cen- 
tralities except the most peripheral bin. This behavior 
might be related to nonflow since many known nonflow 
effects lead to correlations preferentially between oppo- 
site sign particles; e.g. neutral resonances decay into op- 
posite sign particles and jet fragments tend to be charge 
ordered [55]. The LS results therefore typically contain 
smaller nonflow correlations. Bose-Einstein correlations 
between identical particles, on the other hand, can lead to 
larger nonflow for LS than for CI. Figure [I] (right) shows 
the CI and LS results for Cu+Cu collisions at ^/snn =200 
and 62.4 GeV. The same trends hold with the LS results 
lower than the CI results. Figure [2] shows the difference of 
CI t>2{2} and LS «2{2} for Au+Au and Cu+Cu collisions 
at 200 and 62.4 GeV. This difference shows a dependence 
on energy only for central Au+Au collisions. In the low- 
est multiplicity data, CI V2{2} becomes smaller than LS 
v 2 {2}, consistent with expectations from Bose-Einstein 
correlations. 

Figure shows the four-particle cumulant ^2{4} 4 for 
Au+Au (left) and Cu+Cu (right) collisions at 200 and 
62.4 GeV. In the case of v 2 {4:} 4 , no differences are de- 
tected between LS and CI results. This suggests that 
nonflow correlations are suppressed as expected in the 
four-particle cumulant results. Any nonflow source lead- 
ing to fewer than four correlated particles will not con- 
tribute to i>2{4} 4 . In addition, while any nonflow for 
v 2 {2} 2 is suppressed only by 1/M, any nonflow corre- 
lations between four or more particles will still be sup- 
pressed by a combinatorial factor of (M — 1)(M — 2)(M— 
3). «2{4} 4 shows slightly negative values for the more 
central events for Au+Au and Cu+Cu collisions at 200 
and 62.4 GeV. t>2{4} 4 is allowed to take on negative val- 
ues. These may be associated with v 2 fluctuations larger 
than those expected from eccentricity fluctuations alone. 
In this case however, the second or fourth roots of t>2{4} 4 
cannot be defined. For this reason, those points are not 
included in the analysis of v 2 {2} 2 — W2{4} 2 . All results 
are reported in the data tables [53]. It had been ob- 
served from simulations that the measurement of i>2{4} 
using the Q-Cumulants method deviates from input for 
the most peripheral collisions. Also, the LS ^{4} data 
appears to scatter for mean charged particle multiplicity 
density (dN^/dri) < 26. Therefore, no data points are 
used for comparison with models for (dN^/dr]) < 26. 

Figure [I] shows v 2 {2} 2 -v 2 {4} 2 for Au+Au and Cu+Cu 
collisions at 200 GeV (left) and 62.4 GeV (right) for both 
LS and CI. The difference between v 2 {2} 2 and v 2 {&} 2 
is of interest because it is related to nonflow S 2 and v 2 
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FIG. 1: Left: The two-particle cumulant V2{2} 2 for Au+Au collisions at 200 and 62.4 GeV. Results are shown with like-sign 
combinations (LS) and charge-independent results (CI) for 0.15 < pr < 2.0 GeV/c. Right: The same as the left but for Cu+Cu 
collisions. The systematic errors are shown as thin lines with wide caps at the ends and statistical errors are shown as thick 
lines with small caps at the end. Statistical and systematic errors are very small. 
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FIG. 2: The difference of charge-independent (CI) t>2{2} and 
like-sign (LS) i>2{2} for Au+Au and Cu+Cu collisions at 200 
and 62.4 GeV vs. the log of {dN^/drj). The statistical errors 
are smaller than the marker size and not visible for most of 
the data. 



fluctuations: 



v 2 {2} 2 - v 2 {AY ~S 2 + 2< = o 



(8) 



This difference can be taken as an approximate upper 
limit on nonfiow d 2 . We estimate that the approxima- 
tion in Eq. [8] which assumes (v 2 ) is much larger than 
the second, third and fourth moments of v 2 is accurate 
to within 30% for these data sets. We arrive at this 
estimate by assuming v 2 oc e par t and then using our 
Monte Carlo Glauber model to calculate (e pa rt{2} 2 — 
£part{4} 2 )/2<7f part . If the approximation in Eq. [8]is accu- 
rate, this ratio should be unity. We find that for the cen- 
tralizes considered here, the ratio is within 30% of unity. 
Below, where we compare our data to eccentricity mod- 



els, a significant fraction should cancel since the approxi- 
mation applies to both the data and the models. The dif- 
ference v 2 {2} 2 — v 2 {4} 2 increases with beam energy and 
decreases with increasing mean multiplicity. The con- 
tribution from nonfiow typically scales as l/{dN c h,/drj) 
if the number of clusters scales with {dN^/ drj) and the 
number of particles per cluster is constant. A l/JV part de- 
pendence is also expected for a 2 2 from eccentricity fluc- 
tuations. The energy dependence can come from either 
an increase in nonfiow correlations with energy and/or an 
increase in v 2 fluctuations with energy. The LS results 
are systematically lower than the CI results for all but 
the lowest multiplicities, consistent with a nonfiow con- 
tribution to the CI «2{2} results which is reduced for the 
LS i>2{2} results. In the model comparisons that follow, 
we will use the LS results to compare our results to three 
eccentricity models. 



V. DATA AND ECCENTRICITY MODELS 

We compare our v 2 {2} and v 2 {4} results characteriz- 
ing the distribution of v 2 , to equivalent measures char- 
acterizing the eccentricity distributions of three models. 
The models are a Monte-Carlo Glauber model with nu- 
cleons as participants (MCG-N), a Monte-Carlo Glauber 
model with quarks as participants (MCG-Q), and a CGC 
based Monte-Carlo model (fKLN-CGC). The models are 
described in more detail in Appendix [B] Another analy- 
sis of models has been published in Rcf . [30 . The non- 
sphericity of the Au nuclei has been neglected in eccen- 
tricity calculations for the models because non-sphericity 
only affects the most central collisions which are not used 
in the comparison of data with models |31j . 
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FIG. 3: Left: The LS and CI four-particle cumulant v 2 {4} 4 for Au+Au collisions at 200 and 62.4 GeV for 0.15 < p T < 
2.0 GeV/c. The systematic errors are shown as narrow lines with wide caps at the end and statistical errors are shown as thick 
lines with narrow caps at the end. Statistical errors are not visible for most of the points. Right: The LS and CI four-particle 
cumulant t>2{4} 4 for Cu+Cu collisions at 200 and 62.4 GeV for 0.15 < pr < 2.0 GeV/c. The most central points (two points 
for Cu+Cu 62.4 GeV) gives «2{4} 4 < for all the data sets. The negative values are probably due to large fluctuations in 
agreement with Eq. (1). These may include contributions from impact parameter spread and finite multiplicity bin width. 
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FIG. 4: Left: The difference between v 2 {2} 2 and « 2 {4} 2 for 200 GeV Au+Au and Cu+Cu collisions for both LS and CI 
combinations. Right: The difference between V2{2} 2 and ^2{4} 2 for 62.4 GeV Au+Au and Cu+Cu collisions for both LS and 
CI combinations. The statistical and systematic errors are shown as in previous figures. 



A. Upper Limit on Relative Fluctuations 



and contradictory to studies of the nonflow effect. The 
quantity 



We would like to compare our data to models for 
eccentricity fluctuations by comparing <r V2 /v2 to a e je. 
We can not uniquely determine the value of a V2 from 
the two- and four-particle cumulant data however, since 
t>2{2} 2 — «2{4} 2 ~ 82 + 2cr 2 2 . We can however derive an 
upper limit on the ratio cr V2 /v2 by setting 8% = 0. This 
amounts to assuming the difference between the two- and 
four-particle cumulant is dominated by «2 fluctuations 
and that 82 cannot be negative. Although negative non- 
flow values can easily be generated from resonance decays 
in specific kinematic regions, we consider the case that 
the total nonflow should become negative highly unlikely 



R 



«(2-4) 



^ 2 {2} 2 -^ 2 {4p 

V2{2} 2 +V2^Y 



(9) 



then becomes an upper limit to the ratio a V2 / '(^2) where, 
in the case that V2 fluctuations are dominated by eccen- 
tricity fluctuations, (1*2) is the average i>2 relative to the 
participant axis |17) . Additional fluctuations from an- 
other source will lead to a contribution to the difference 
between i>2{2} and ^{4} not related to the eccentricity 
fluctuations that relate the reaction plane and the par- 
ticipant plane. In the following figures, we compare the 



ratio i?„(2-4) for the like-sign results to the ratio 



R 



■e(2-4) 



/e{2} 2 -e{4} 2 
e{2} 2 + e{4} 2 



(10) 



for the three eccentricity models described in Ap- 
pendixjBj where e{2} and e{4} are the second and fourth 
cumulants for e par t . Since higher moments (skewness and 
kurtosis) of the distribution of V2 or £2 contribute to Eqs. 
[9] and [lO] it is important to compare the same quantities 
from data and the eccentricity models. For small e Eq.[T0| 
becomes a e /(e). 

Figure [5| shows R v (2-4) vs mean charged hadron mul- 
tiplicity for 200 GeV (left) and 62.4 GeV (right) Au+Au 
data. The LS i<2{2} results are used to reduce nonflow. 
The data is compared to the same quantity for the three 
different models. The shaded bands show the uncertain- 
ties on the models that arise primarily from the uncer- 
tainty in the Woods-Saxon parameters used to describe 
the nuclei. The error is correlated between Monte-Carlo 
models and for clarity is only plotted on the MCG-Q 
and fKLN-CGC models. The centrality in the models is 
defined using multiplicity so that the model calculations 
include bin-width effects and impact parameter fluctua- 
tions similar to data. In as much as the models correctly 
model the multiplicity, by defining centrality in the mod- 
els the same way that it is defined in data, both the model 
and the data will have the same impact parameter fluc- 
tuations. 

In peripheral collisions ((dN c h/ drj) < 150), data ex- 
ceeds the eccentricity models substantially. This is not 
surprising since we expect a significant contribution from 
nonflow in this region. The central value for the ratio 
from the MCG-N model rises with increasing centrality 
and then overshoots the upper limit in the most central 
collisions. Given the errors indicated by the yellow band 
however, the MCG-N model could still be consistent with 
the upper limit. The MCG-Q model approaches the up- 
per limit in central collisions but never exceeds it. The 
fKLN-CGC model has the smallest values and is well 
below the upper limit throughout the entire centrality 
range. Notice that in the models, the more constituents, 
the smaller the fluctuations. 

In Fig. [5] (right), the 62.4 GeV Au+Au data are com- 
pared to models. Data points are only reported where 
«2{4} 4 is positive. At this lower energy, peripheral data 
is again above the models. The central value for the 
MCG-N model again overshoots the upper limit for cen- 
tral and mid-central collisions while the MCG-Q model 
appears to just reach the upper limit for the most central 
data point. The uncertainty on the geometry of the Au 
nucleus again however makes it impossible to rule out any 
of the models in this comparison. The fKLN-CGC model 
lies below the upper limit for the entire range. The fact 
that the MCG-N and MCG-Q models reach and in some 
cases exceed the upper limit means that for those models 
to be correct, nonflow would have to be small or perhaps 
even negative. Nonflow can be negative from resonance 



decay, but is not likely. The lower energy data therefore 
provides a very useful test of the models and results from 
the beam energy scan at RHIC promise to provide even 
better constraints |32j . 

Figure [6] shows the STAR 200 GeV Au+Au data on 
the upper limit for cr V2 /(v2} compared to the PHOBOS 
results reported in Ref. [33 under their assumption that 
62 is zero for Arj > 2 (see the reference for details). The 
PHOBOS results are for all charged particles while the 
STAR results are for LS pairs only. PHOBOS has sub- 
tracted narrow A77 correlations by fitting W2(^i)w2(%) 
and removing the narrow diagonal peak corresponding 
to small- A77 nonflow correlations. This may explain why 
the PHOBOS results are slightly below the STAR upper 
limits derived from LS «2{2}, suggesting that there may 
be some residual nonflow in our LS results. We also note 
however, that the analysis procedures in this paper and 
in the PHOBOS paper are quite different. 

Figure [7] shows the upper limits and models for Cu+Cu 
collisions at 200 and 62.4 GeV (respectively left and 
right). Data points are only reported where i>2{4} 4 is 
positive. The upper limit on fluctuations for Cu+Cu col- 
lisions are larger than for Au+Au and lie near unity. All 
the models fall below the upper limit and differences be- 
tween the models are small. This is likely due to the large 
multiplicity fluctuations for smaller systems in the mod- 
els which masks the other physical differences between 
the models. The large Cu+Cu results do not provide 
constraint on the models. 



B. Nonflow 

An alternative way to look at the t>2{2} 2 and i>2{4} 2 
data is to apply the assumption that V2 fluctuations are 
dominated by the initial spatial eccentricity fluctuations 



"2 



(11) 



and then derive the nonflow 82 that would be implied by 
each eccentricity model, 
observable. 



In Eq. 



11 



«2 



is not directly 

Then we can calculate the value of 62 that 



would be needed to satisfy the following equation: 
v 2 {2} 2 -v 2 {4} 2 ~6 2 + 2a 2 V2 . 



(12) 



Recalling from Eq. [T] and [2] that v 2 {2} 2 + u 2 {4} 2 S 2 
2(«2) 2 , we derive the following expression for 82' 



82 « v 2 {2} 2 - v 2 {4} 2 



(13) 



which only depends on the directly observable cumulants 
and quantities obtained from models. Since a model de- 
pendence exists, the 82 values are not measurements of 
82 but instead provide an alternative consistency check 
for the models. These values can be compared to other 
measurements of nonflow correlations such as the already 
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FIG. 5: The upper limit on a V2 / ' (V2) for 200 GeV (left) and 62.4 GeV (right) Au+Au collisions from Eq. (|9| compared to a e /e 
from Eq. (10 1 for three different models. The upper limit is found using the LS results for «2{2}. Data are from the range 
The shaded bands reflect the uncertainties on the models which are dominated by uncertainty on 
The uncertainty is only shown for the MCG-N and fKLN-CGC models. The 



0.15 < p T < 2.0 GeV/c 

the distribution of nucleons inside the nucleus 
uncertainty on the MCG-Q model is the same as for the MCG-N model but is not shown for the visual clarity. 
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FIG. 6: The STAR data compared to PHOBOS data [33] on 
&v 2 /{ v 2) with 82 for Arj > 2 taken to be zero (see Fig. 6 from 
Ref. 

Ref. [33]. 



]). The shaded band shows the errors quoted from 



measured two-particle correlations [21j . This is an impor- 
tant test for the models, since a complete model of heavy 
ion collisions should be able to predict multiple observ- 
ables at once. The interpretation, however, of the struc- 
tures in two-particle correlations such as the ridge [21] is 
in flux. In particular, the nonflow correlations from jets 
are inferred from two-particle correlations vs Ar) and A(j) 
after subtracting a Ar/ independent v\ term. This ap- 
proximation may not be valid for reasons discussed re- 
cently in the literature [34]. Given the current state of 
understanding, in this paper we do not make a direct 
comparison of the nonflow correlations inferred from this 
analysis to those inferred from two-particle correlations. 



In the absence of new physics, the term 82 will vary 
with event multiplicity as 1/M. This is because, in the 
case that high multiplicity events are a linear superpo- 
sition of lower multiplicity events, the numerator in the 
mean grows as M while the denominator grows as the 
number of pairs M(M — l)/2. To cancel out the combi- 
natorial 1/M dependence we scale 82 by the number of 
mean charged hadrons within < 0.5. A variation of 
( rf ^ h )c?2 with multiplicity implies a non-trivial change in 
the physics. 

Figure [8j (left) shows the like-sign (^ k )8 2 that is re- 
quired if the Monte-Carlo Glauber model with nucleon 
participants gives the correct description of the eccentric- 
ity fluctuations and eccentricity fluctuations dominate i>2 
fluctuations. The nonflow is larger at 200 GeV than at 
62.4 GeV. Within errors, (^)<5 2 is the same in Cu+Cu 
collisions and Au+Au collisions at the same energies and 
event multiplicities. The value of ( C ^/ L )82 required by 
this model seems to fall with centrality and the central 
value becomes negative for the most central Au+Au col- 
lisions. The errors shown in the figure are dominated 
by the systematic errors on the MCG-N model which 
are highly correlated from point to point (they depend 
on the parameters for the geometric description of a Au 
nucleus). As such, while the most central data point is 
still consistent with zero, the dropping trend with cen- 
trality is significant for most of the range allowed for 
describing Au nucleus. For this model of eccentricity 
fluctuations to be valid, the nonflow in central Au+Au 
collisions would have to be near zero or negative. This 
appears to contradict measurements showing significant 
near-side two-particle correlations [21] [35] . In case other 
sources besides eccentricity fluctuations contribute to V2 
fluctuations, the inferred nonflow would need to become 
even smaller. While the uncertainties arising from the 
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FIG. 7: The upper limit on a V2 /(v 2 ) for 200 GeV (left) and 62.4 GeV (right) Cu+Cu collisions from Eq. |9| compared to a e /e 
from Eq. ( 10 1 for three different models. 
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FIG. 8: The mean multiplicity scaled nonflow (62) found by assuming Eq. 11 and using the three different models for a 2 /e 2 : 
MCG-N (left), MCG-Q (middle), and fKLN-CGC (right). The systematic errors are shown as thin lines with wide caps and 
statistical errors are shown as thick lines with narrow caps. Statistical errors are not visible for most of the points. 



geometric description of the Au nucleus preclude a defini- 
tive statement, it appears likely that the MCG-N model 
over-predicts the ratio of eccentricity fluctuations to the 
mean eccentricity. 

Figure 8 (middle) shows the {^^ L )8 2 required if the 
Monte-Carlo Glauber model with constituent quark par- 
ticipants gives the correct description of the eccentric- 
ity fluctuations and if eccentricity fluctuations dominate 
v 2 fluctuations. Within errors, (^g— )5 2 is the same in 
Cu+Cu collisions and Au+Au collisions at the same ener- 
gies and event multiplicities. The smaller relative fluctu- 
ations for the constituent quark participant model means 
this model would be consistent with larger nonflow values 
than the nucleon participant model. The required non- 
flow values are essentially positive at all measured mul- 
tiplicities. This means this model has a better chance 
of accommodating the near-side two-particle correlations 
observed in data. 

Figure 8 (right) shows (^^ k )5 2 derived using the 
fKLN-CGU Monte-Carlo model. This model has a larger 
average eccentricity and smaller eccentricity fluctuations 



leading to the smallest relative fluctuations of the three 
models. The mean multiplicity scaled nonflow again is 
larger for 200 GeV collisions than 62.4 GeV collisions and 
Cu+Cu collisions seem to have the same nonflow values 
as Au+Au when they are compared at the same mean 
multiplicity. The multiplicity scaled nonflow implied by 
the fKLN-CGC eccentricity model increases slightly or 
remains flat with centrality. CGC models for the ini- 
tial conditions of heavy-ion collisions have also been in- 
voked to try to explain the near-side correlations ob- 
served in data [36| . This analysis adds information from 
four-particle correlations not accessible through measure- 
ments of a two-particle correlation function. It remains 
to be seen if a consistent determination of two- and four- 
particle cumulants related to v 2 , v 2 fluctuations and non- 
flow can be derived from a CGC model with radially 
boosted flux tubes. 
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C. Eccentricity Scaling of v% 

Now we show the ratio {v2)/{s) for the three models 
of eccentricity. In the case that v-i oc e, then (v2)/(s) in 
the reaction plane is given by i; 2 {4}/e{4} [3"7] (e{4} is 
the fourth cumulant defined in Appendix [b]) . In Fig. [9] 
we plot W2{4}/e{4} vs mean multiplicity forAu+Au and 
Cu+Cu collisions at 200 and 62.4 GeV. When plotted 
vs mean multiplicity, v-i oc e from all systems and ener- 
gies seem to fall on top of each other. The fKLN-CGC 
model displays a saturation with V2 cc e. The Monte 
Carlo Glauber model with Nucleon participants shows 
the steepest increase of vi/e while the constituent quark 
model is intermediate between the sharp rise of the nu- 
cleon participant model and the saturation of the fKLN- 
CGC model. The approximation that t>2 oc e is strongly 
violated for the nucleon participant model. This also im- 
plies that ^2{4}/e{4} = (V2) / (s) may be broken since 
that equality holds only when vi oc e. The violation of 
V2 oc e also implies that if the nucleon participant model 
is the correct eccentricity model, then the collisions at 
RHIC may be far from the ideal hydrodynamic limit. 
The fKLN-CGC model and constituent quark model im- 
ply x>2 saturates or nearly saturates in central Au+Au 
collisions, consistent with a nearly perfect liquid behav- 
ior. 



correlations data. We presented V2/e for the three differ- 
ent eccentricity models and found that the fKLN-CGC 
model for eccentricity leads to a saturation of v-ile for 
Au+Au collisions with (^^+) > 300 while v 2 /e is ris- 
ing at all centralities when the MCG-N model is used for 
e. Assuming fKLN-CGC to describe the initial state ec- 
centricity, the saturation of V2/S provides support for a 
nearly perfect hydrodynamic behavior for heavy-ion col- 
lisions at RHIC. 
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VI. CONCLUSIONS 



Appendix A: Q-Cumulants vs Fitting q-distributions 



We presented STAR measurements of two- and four- 
particle V2 cumulants (t>2{2} and t>2{4}) for Au+Au and 
Cu+Cu collisions at ^/s NN = 200 and 62.4 GeV along 
with the difference v 2 {2} 2 — v 2 {<i} 2 5 2 + 2a V2 = a 2 ot 
for charge-independent and like-sign combinations of par- 
ticles. «2{4} 4 shows negative values for the most cen- 
tral collisions for all the data sets. The difference 
w 2 {2} 2 — w 2 {4} 2 increases with beam energy for both 
Cu+Cu and Au+Au collisions. For a given ^/snn and 
mean charged particle multiplicity, t>2{2} 2 — «2{4} 2 val- 
ues are the same in Cu+Cu and Au+Au collisions within 
errors. Although the value of i>2 fluctuations can not 
be uniquely determined in this way, f2{2} and t>2{4} 
were used to place an upper limit on the ratio cr„ 2 /w 2 . 
The eccentricity fluctuations from the MCG-N model are 
largest, rising above the upper limit from data for cen- 
tral Au+Au collisions, but the MCG-Q and fKLN-CGC 
eccentricity models fall within the presented limit. To 
further investigate the models we calculated the value 
of the nonflow 62 implied by the models for eccentricity 
fluctuations under the assumption that cr V2 /v2 = cr e /e. 
The nonflow values implied by the fKLN-CGC model 
are larger than those from either of the Monte Carlo 
Glauber models. The nonflow implied by the fluctuations 
in the MCG models becomes zero or negative for central 
Au+Au collisions. This analysis challenges theoretical 
models of heavy-ion collisions to describe all features of 
the data including V2, v 2 fluctuations, and the various 



The fitting of the reduced flow vector distribution is 
described in more detail in Ref. [TH]. The fit param- 
eters described in that reference can be transformed to 
w 2 {2,qfit} 2 ee U2 {qfit} 2 +c7 2 ot and u 2 {4,qfit} 2 = v 2 {q&t} 2 
where «2{2,qfit} and t>2{4, qfit} are the two- and four- 
particle cumulants determined from the q-distribution 
which can be compared to other determinations of v 2 {2} 
and f2{4}. In Fig. 10 (top) we show the ratio of «2{2} 
determined from the q-distribution analysis and the Q- 
Cumulants analysis. 

Deviations between the q-distribution and Q- 
Cumulants results can be seen when the multiplicity of 
the event is smaller with the q-distribution results being 
smaller than the Q-Cumulants results. These deviations 
can be traced to the break-down of the large N approx- 
imation required when fitting the q-distribution. An 
attempt is made to correct for this break-down which 
brings the results closer together but the deviations 
are still significant for multiplicities below 150. The 
correction is carried out by adjusting the q-distribution 
data before it is fit. The correction is derived by taking 
the ratio of the expected and observed q-distribution 
from simulated data. Although the correction extends 
the apparent validity of the q-distribution analysis to 
lower multiplicities, we find that the q-distribution 
analysis is less reliable than the Q-Cumulants analysis. 



Figure 10 (bottom) shows the ratio of the quantity 
«2{2} 2 — wTH} 2 from the q-distribution fits over the same 
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■A- Cu+Cu 62.4 GeV 
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FIG. 9: The eccentricity scaled V2 for 200 and 62.4 GeV Au+Au and Cu+Cu collisions with eccentricity taken from the MCG-N 
(left), MCG-Q (middle), or fKLN-CGC (right) model. The statistical and systematic errors are shown as in previous figures. 
Statistical errors are not visible for most of the points. 




0.6 



100 



200 



300 

<dN /dn> 

cn 



400 



500 



600 



FIG. 10: Top panel: the ratio of the two-particle cumulant 
U2{2} for Au+Au collisions at 200 and 62.4 GeV evaluated us- 
ing the q-distribution method and the Q-Cumulants method. 
Both results are calculated for combinations of particles in- 
dependent of their charge (CI). Bottom panel: The ratio of 
the q-distribution and the Q-Cumulants method results for 
v 2 {2} 2 - v 2 {4} 2 (CI) for 200 GeV Au+Au and Cu+Cu col- 
lisions. In both panels, systematic errors are shown as thin 
lines with wide caps and statistical errors are shown as thick 
lines with narrow caps. 



from the Q-Cumulants analysis. Data are from 200 GeV 
Au+Au and Cu+Cu collisions. The two methods pro- 
duce significantly different results for i> 2 {2} 2 — «2{4} 2 
with the difference most pronounced in Cu+Cu and 



peripheral Au+Au collisions. The q-distribution gives 
smaller values. This is related to the large N approx- 
imation required in the fitting procedure for the q- 
distribution. When multiplicity is low, the tails of the 
q-distribution cannot be populated. We find that this 
leads to a narrowing of the observed distribution relative 
to the fit function and the width of the distribution deter- 
mines V2{2} 2 — «2{4} 2 . The q-distribution fits therefore 
underestimate V2{2} 2 — 1>2{4} 2 , so we use the results from 
the Q-Cumulants calculation in this paper. 



Appendix B: Three Eccentricity Models 

We use three Monte-Carlo models to study eccentricity 
and eccentricity fluctuations. The first two are Glauber 
models which either treat nucleons as participants or 
constituent quarks within the nucleons as participants 
(MCG-N and MCG-Q respectively). The third model is 
the factorized Kharzeev, Levin, and Nardi Color Glass 
Condensate model (fKLN-CGC) [131. The input param- 
eters used for the Woods-Saxon distribution of nucleons 
are in Table II. The Au nuclei have been assumed spher- 
ical for the eccentricity calculations. A 0.4 fm exclusion 
radius is used in the calculations so that nucleons do not 
overlap in coordinate space. The MCG-N model is de- 
scribed elsewhere [3 [38] and is used to calculate the iVp ar t 
and N hin values in Table [TIT} For the MCG-Q model, 
we first distribute nucleons inside a nucleus according 
to a Woods-Saxon distribution with parameters taken 
from [39] , then we distribute three constituent quarks in- 
side each nucleon according to another Woods-Saxon dis- 
tribution where the radius of the nucleon is taken to be 
0.63 fm and the surface width is 0.08 fm. The results were 
not very sensitive to variations of these parameters within 
a reasonable range. One might consider a Gaussian for 
the quarks instead of a Woods-Saxon distribution. The 
Woods-Saxon distribution gives a more flat-topped dis- 
tribution but the calculated eccentricity and eccentricity 
fluctuations are not highly sensitive to the exact distri- 
bution. The main feature of the MCG-Q model is that 
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the potential number of participants increases by a factor 
of three and there are large correlations between partici- 
pants because the quarks are confined within the nucle- 
ons. 



Parameter / System 


197 Au + i!)Y Au 63 Cu + S3 Cu 


R 

a 


6.38 ±0.06 fm 4.218 ± 0.014 fm 
0.535 ± 0.027 fm 0.596 ± 0.005 fm 



TABLE II: Input parameters for Woods-Saxon distribution in 
Monte-Carlo Models. 

The Woods-Saxon parameters from f3"9l are based on 
measurements of electron scattering which are only sen- 
sitive to protons. If the Au nucleus has a neutron skin, 
then the hadronic radius may be larger than that quoted 
in [3j3]. We estimated the systematic errors by varying 
the Woods-Saxon parameters within the range allowed by 
electron scattering data. Although unmeasured, theoret- 
ical guidance suggests the neutron skin may add 0.2 fm to 
the radius of heavy nuclei j3Hj ■ To account for a possible 
neutron skin, we increased the radius of the Au nucleus 
to 6.7 fm. We find that our results only weakly depend 
on the radius and depend mostly on the diffuseness pa- 
rameter " a" . The effect of a neutron skin is therefore well 
within our quoted systematic errors and will not affect 
our conclusions unless the skin significantly changes the 
diffuseness at the edge of the nucleus. 

The fKLN-CGC model provides multiplicity and ec- 
centricity. Our MCG calculations use a two-component 
model and a negative binomial distribution to estimate 
the event multiplicity for each simulated event. The first 
parameter of the binomial distribution is generated for 
each event using 



3+ard 

) + 2x h ardA r bin/-/Vpart) (Bl) 

0.5933 ln(ys^/GeV/c 2 ) - 0.4153, 




part 



where /( 

A^bin is the number of nucleon-nucleon collisions, 
is the number of participating nucleons, and Xhard is t ne 
fraction of the multiplicity proportional to N^- m . Then 
multiplicity is generated by sampling a negative bino- 
mial distribution with parameters n and width k = 2.1 
for each participant. This parametrization provides a 
good description of multiplicity measurements in heavy- 
ion collisions from ^/s NN = 20 to 200 GeV [H] and for 
all centralities. For the MCG-Q model, while the eccen- 
tricity is defined by the locations of constituent quarks 
participating in the collisions, the multiplicity is defined 
by the nucleon Ap art and A^in- We define the central- 
ity of the models according to this multiplicity so that 
the data and model are treated equivalently. In this way, 
our eccentricity fluctuations also contain the impact pa- 
rameter and A^part fluctuations that are intrinsic to our 
experimental determination of a given centrality interval. 
The uncertainties on the models were estimated by vary- 
ing the Woods-Saxon parameters within the range of the 
errors quoted in ref. |39j . We also varied the parameters 



for the multiplicity but the results were not very sensitive 
to those. 
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FIG. 11: e st d (top) and a e (bottom) vs. centrality for Au+Au 
200 GeV between Monte Carlo Glauber-Nucleon Participants, 
Monte Carlo Glauber-Quark Constituents, Color Glass Con- 
densate models. The shaded regions show the systematic er- 
rors. 

Several different variables related to the eccentric- 
ity have been calculated from the three models. This 
includes the eccentricity relative to the reaction-plane 
(^std = U)+x) )i the- eccentricity relative to the participant 



^std 

plane (e part 
of e 



and the two- and four-particle cumulants 



part 



e{2} = 

£ {4} = (2<4 art 



part / 
2 \2 



- part 



)) 



1/4 



(B2) 
(B3) 



e s td for 200 GeV Au+Au collisions is shown in Fig. [TT] 
(top). e st d is largest for the fKLN-CGC model and small- 
est in the MCG-N model. The MCG-Q model is inter- 
mediate between the two. The relevant quantities have 
been tabulated online |29j . 

Figure 11 (bottom) shows the fluctuations of £ s td for 
the three models for 200 GeV Au+Au collisions. The 
fluctuations in the two Glauber models are larger than 
those for the fKLN-CGC model. One might expect the 
MCG-Q model to have smaller fluctuations than the 
MCG-N model since there are three times as many possi- 
ble participants. This is counterbalanced however by two 
effects 1) the three constituent quarks are confined inside 
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nucleons thus inducing correlations that partially offset 
the effect of more participants, and 2) the mean value 
of the eccentricity is larger in the MCG-Q model. These 
effects lead to the result that the width of the eccentric- 
ity distribution in the MCG-Q model and the MCG-N 
model are similar. On the other hand, since the MCG-Q 
model gives a larger average eccentricity, when consider- 
ing cr e /e, the MCG-Q model is intermediate between the 
fKLN-CGC and the MCG-N models as one might naively 
expect. 

The trends for Cu+Cu collisions remain the same as 
for Au+Au collisions with the fKLN-CGC model having 
the largest eccentricity and smallest fluctuations and the 
MCG-Q model intermediate between the MCG-N and 
fKLN-CGC models. None of the models showed a signif- 
icant difference between A /s NN = 62.4 and 200 GeV, so 
we only show the 200 GeV results here. 



ripheral collisions. MCG-Q and fKLN-CGC models show 
the same values for e{2} 2 for central and mid-central col- 
lisions but MCG-Q shows the highest values in all the 
three models for the most peripheral collisions. e{4} 4 
shows similar behavior as e{2} 2 but it becomes negative 
for the most central collisions in all the models like t>2{4} 4 
for the most central collisions in the data. This behavior 
is the same for Cu+Cu collisions and different energies. 
In the models, this negative value can be traced to A part 
fluctuations present when using multiplicity to select cen- 
trality bins. If A part is used to define the centrality in 
the models, then e{4} 4 remains positive, even for central 
collisions. 

TABLE III: The (dN^/drj) [12], iVp art and 7V bin values cor- 
responding to the centrality intervals used in this paper. 
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FIG. 12: e{2} 2 and e{4} 4 



centrality for Monte Carlo 



Glauber models with nucleon or constituent quark partici- 
pants and for a Color Glass Condensate model. 



Au 



Figure 12 shows e{2} (top) and e{4} (bottom) for 
Au 200 GeV for the three models. e{2} 2 shows pos- 
itive values throughout the range and decreases with in- 
creasing centrality. The MCG-N model shows smaller 
values than the other two models for central and mid- 
central collisions but cross fKLN-CGC for the most pe- 
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Au+Au 200 GeV 



22 ±2 
45 ±3 
78 ±6 
126 ±9 
195 ±14 
287 ±20 
421 ±30 
558 ±40 
691 ±49 



13.46 ±0.50 
26.62± 0.95 
47.06 ±1.21 
75.58±1.56 
114.81±1.73 
166.85±1.33 
234.49±0.84 
299.47±0.75 
349.09±0.30 



12.45 ± 0.69 
29.33 ± 1.75 
62.1 ± 2.1 
121. 8± 4.2 
218.9 ± 6.1 
371.3 ± 6.2 
599.6 ± 4.5 
845.6 ± 3.2 
1059 ± 3. 



Au+Au 

13.9 ±1.1 
29.1 ±2.2 

53.1 ±4.2 

87.2 ±7.1 
135 ±11 
202 ±17 
292 ±25 
385 ±33 
472 ±41 
Cu+Cu 



62.4 GeV 

13.18 ±0.71 
25.56± 1.11 
44.97 ±1.27 
72.70±1.25 
110.53±1.05 
161.08±0.97 
228.51±0.52 
293.39±0.96 
343.82±0.44 
200 GeV 



11.6 ± 0.87 
26.69 ± 1.87 
55.9± 2.9 



107.4 
190.5 
319.4 
514.6 
721.7 
900.75 



3.8 
4.5 
4.8 
3.4 
3.9 
1.85 



25.3 ±1.6 
38.7 ±2.5 
56.9 ±3.7 
82.9 ±5.4 
119 ±7.7 
170 ±11 



16.41±0.24 
25.14±0.16 
37.35±0.47 
53.07±0.29 
73.61±0.12 
98.08±0.11 



15.71±0.31 
27.42 ± 0.22 
46.87 ± 1.00 
75.46 ± 0.42 
119.65 ±0.15 
182.7 ± 0.30 



Cu+Cu 62.4 GeV 



17.4 ±1.1 

26.3 ±1.7 
38.7 ±2.5 

56.4 ±3.7 
81.2 ±5.3 
117 ±7.7 



15.36±0.05 
23.92±0.05 
35.62±0.05 
50.76±0.12 
70.67±0.50 
94.98±0.25 



13.91±0.01 
25.56±0.06 
41.09± 0.12 
65.86±0.30 
103.15±0.95 
155.65±0.75 
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